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The PHENIX collaboration has measured (f meson production in d+Au collisions at ,/sffff = 

200 GeV using the dimuon and dielectron decay channels. The (j) meson is measured in the forward 
(backward) d-going (Au-going) direction, 1.2 < y < 2.2 (—2.2 < y < —1.2) in the transverse- 
momentum (pt) range from 1-7 GeV/c, and at midrapidity \y\ < 0.35 in the pr range below 
7 GeV/c. The cji meson invariant yields and nuclear-modification factors as a function of pr, rapidity, 
and centrality are reported. An enhancement of <j) meson production is observed in the Au-going 
direction, while suppression is seen in the d-going direction, and no modification is observed at 
midrapidity relative to the yield in p-\-p collisions scaled by the number of binary collisions. Similar 
behavior was previously observed for inclusive charged hadrons and open heavy flavor indicating 
similar cold-nuclear-matter effects. 

PACS numbers: 25.75.Dw 


I. INTRODUCTION 

Collisions of deuterons with gold nuclei (d-|-Au) are of 
significant interest in the study of the strongly coupled 
Quark Gluon Plasma (QGP) produced at the Relativis¬ 
tic Heavy Ion Collider (RHIC) UM- The highly asym¬ 
metric collisions of a small projectile and a large target 
nucleus provide a way to investigate the initial state of a 
nucleus-nucleus collision experimentally, potentially dis¬ 
entangling effects due to QGP formation from the cold- 
nuclear-matter effects. The latter include modification of 
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the parton distribution functions (PDFs) in the nucleus 
relative to those in the nucleon Q, initial-state energy 
loss 0, and the so-called Cronin effect. The Cronin ef¬ 
fect refers to the enhancement of high-p^ particle produc¬ 
tion in p-\-A collisions relative to that in p-\-p collisions 
scaled by the number of binary collisions and is often 
attributed to multiple scattering of the incoming parton 
inside the target nucleus [1-Q. In addition, results from 
recent p{d)-\-A collisions at the Large Hadron Collider 
and RHIC suggest that long-range correlations, either 
present in the initial state or induced by the evolution of 
the medium, play an important role even in these small 
collision systems [^i3- Detailed studies of particle pro¬ 
duction systematics in c?-|-Au at RHIC may inform this 
question pll - l^ . 

PHENIX has measured the production of identihed 
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particles, such as tt, K andp, at midrapidity in d+Au and 
Au+Au collisions and has found intriguing similarities 
between the AT/tt and p/tt ratios in peripheral Au+Au 
and central d+Au collisions [2^. However, it was ob¬ 
served that the ratio of spectra in peripheral Au+Au to 
those in central d+Au starts above one at low-pT and 
trends to a constant value of ^ 0.65 for all identified par¬ 
ticles at high pt- One explanation for the \ow-pT rise 
is a relative deficit of midrapidity soft particle yield in 
d+Au collisions compared to Au+Au, which could be 
due to the participant asymmetry in the d+Au collisions 
producing a rapidity shift in the peak of particle produc¬ 
tion 0. Consequently, measuring identified particles in 
different regions of rapidity provides a more complete pic¬ 
ture of d+Au collisions and sheds light on the relationship 
between p+p, d+Au and Au+Au. Measuring identified 
particle production at forward rapidity also enables ac¬ 
cess to the low-a: region where nuclear PDFs (nPDFs) are 
not well known and where one expects parton-saturation 
effects to begin to manifest in modified particle produc¬ 
tion. 

The yield of (j) mesons in high energy heavy ion col¬ 
lisions provides key information about the QGP, as the 
yield is potentially sensitive to medium-induced effects 
such as strangeness enhancement , a phenomenon as¬ 
sociated with soft particles in bulk matter which can be 
accessed through the measurements of (j) meson produc¬ 
tion (28l-l3^. PHENIX measures (j) meson production 
over a wide rapidity range and for various collision sys¬ 
tems such as p+p, d+Au and Au+Au. The production 
of (j) mesons has already been measured in p+P, d+Au, 
Cu+Cu and Au+Au at midrapidity [^ . [^ and in 
p+p at forward and backward rapidities over a wide 
range in pT ■ 

In this paper, the production of (j) mesons is determined 
at forward and backward rapidities via dimuons detected 
in the PHENIX muon spectrometers and at midrapidity 
via dielectrons detected in the PHENIX central arms. 
Measurements of the </> meson nuclear-modification fac¬ 
tor (RdAu) in d+Au collisions versus rapidity and versus 
Pt are presented. The results presented here are based 
on d+Au collisions at ,/s^ = 200 GeV recorded in 2008. 
The luminosity used in </) —>■ i +^analysis corresponds 
to 60 nb“^ which is equivalent to a nucleon-nucleon inte¬ 
grated luminosity of 24 pb“^, while the luminosity used 
in ^ e“''e“ analysis corresponds to 71 nb“^ which is 

equivalent to a nucleon-nucleon integrated luminosity of 
28 pb-i. 


II. EXPERIMENT 

Figure [T] shows a schematic of the PHENIX detector, 
which is described in detail in Ref. [d^. The detectors 
relevant for the analysis of the (j) meson in the dilepton 
channels are the central arms, two muon spectrometers 
and the two beam-beam counters. 

The PHENIX central arms measure particles by using 



18.5 m = 60 ft 


FIG. 1. (color online) A schematic side view of the PHENIX 
detector configuration in the 2008 data-taking period. 


drift chambers (DG) and pad chambers for the tracking 
and momentum measurement of charged particles, ring 
imaging Cerenkov detectors (RICH) for the separation of 
electrons up to the tt meson Cerenkov radiation threshold 
at 5 GeV/c, and an electromagnetic calorimeter (EMCal) 
for the measurement of spatial positions and energies of 
photons and electrons. The EMCal comprises six sectors 
of lead-scintillator calorimeter and two sectors of lead- 
glass calorimeter. 

The muon spectrometers, located at forward and back¬ 
ward rapidity, are shielded by absorbers composed of 
19 cm of copper and 60 cm of iron. Each spectrometer 
comprises the muon tracker (MuTr) immersed in a radial 
magnetic field of integrated bending power of 0.8 T-m, 
and backed by the muon identifier (MuID). The muon 
spectrometers subtend 1.2 < \r]\ < 2.2 and the full az¬ 
imuth. The MuTr comprises three sets of cathode strip 
chambers while the MuID comprises five planes of larocci 
tubes interleaved with steel absorber plates. The momen¬ 
tum resolution, 5p/p, of particles in the analyzed momen¬ 
tum range is about 5%, independent of momentum and 
dominated by multiple scattering. Muon candidates are 
identified by matching tracks reconstructed in the MuTr 
to MuID tracks that penetrate through to the last MuID 
plane. The minimum momentum of a muon to reach the 
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last MuID plane is ^2 GeV/c. 

Beam-beam counters (BBC), comprising two arrays of 
64 Cerenkov counters covering the pseudorapidity range 
3.1 < It^I < 3.9, measure the collision vertex along the 
beam axis (zvtx) with 0.5 cm resolution in addition to 
providing a minimum-bias (MB) trigger. 

III. DATA ANALYSIS 

This section describes the details of the measurement 
of the (j) meson in 0 —>■ and (p —>■ e+e" decay 

channels. 

A. 4 > ^ measurement 

The dimuon data set for this analysis was collected in 
2008 using a MB trigger that requires at least one hit in 
each of the BBC detectors in conjunction with the MuID 
Level-1 dimuon trigger. Also, at least two tracks were 
required to penetrate to the last layer of the MuID. The 
MB trigger measures 88 ± 4% of the total d-l-Au inelastic 
cross section of 2260 mb [2l|. 

A set of cuts was employed to select good muon candi¬ 
dates and improve the signal to background ratio. Event 
selection requires the BBC collision z vertex to be re¬ 
constructed within ±30 cm of the center of the interac¬ 
tion region along the beam axis. The MuTr tracks are 
matched to the MuID tracks at the first MuID layer in 
both position and angle. The track is required to have 
at least 8 of 10 possible hits in the MuID. 

The invariant mass distribution is formed by combin¬ 
ing muon candidate tracks of opposite charge. In ad¬ 
dition to low mass vector mesons, the invariant mass 
spectra in the region of interest contain uncorrelated and 
correlated backgrounds. The uncorrelated backgrounds 
arise from random combinatoric associations of unrelated 
muon candidates while the correlated backgrounds arise 
from open charm decay (e.g., DD where both mesons 
decay semileptonically to muons), open bottom decay, rj 
and ui Dalitz decays and from the Drell-Yan process. The 
correlated background is much smaller than the uncor¬ 
related background for all centralities. It is also estab¬ 
lished from simulation that the background is dominated 
by the uncorrelated contribution from decays, such as 
AT, TT —>■ fiv, that occur in front of and inside the ab¬ 
sorber. A smaller contribution to the background comes 
from hadron decays in the muon tracker volume. When 
combined into track pairs, these pairs produce a broad 
distribution of invariant masses and also a broad distri¬ 
bution of XvtxJ parameter resulting from fitting the 
two muon tracks with the BBC measured event vertex 
position. This procedure separates the foreground and 
background spectra by applying a cut of Xvtx ± 4.0 to 
extract foreground spectra and a cut of Xvtx > 4.0 to ex¬ 
tract background spectra. The cut value, xltx cut = 4.0, 
is selected such that it retains as much of the signal as 


possible within the foreground, while still allowing suffi¬ 
cient statistics in the background sample. This procedure 
is described in detail in Ref. [s^. 

To extract the (p meson signal, the background pairs 
are subtracted from the foreground pairs. Because the 
pairs with Xvtx ± 4.0 represent only the uncorrelated 
part of the background the shape of their mass distri¬ 
bution is slightly altered to account for the contribution 
of correlated pairs, which have much smaller contribu¬ 
tion according to simulation results. This is done by 
fitting a fourth order polynomial to the ratio of the fore¬ 
ground to background distributions in the nonresonance 
region (0.3 < M^+^- < 0.5 GeV/c^ and 1.5 < M^+^- < 
2.5 GeV/c^). The background distribution is then mul¬ 
tiplied by the polynomial in the entire mass range. Fig¬ 
ure [5] shows the resulting foreground and renormalized 
background distributions. 



FIG. 2. (color online) The unlike-sign dimuon invariant mass 
spectrum (solid red points) and the renormalized background 
spectrum (open blue circles). 

In the mass range 0.3 < M^+^- < 2.5 GeV/c^, the 
resulting distribution after subtraction has contributions 
from three mesons: w, p, and p mesons. The distribution 
is fitted with the sum of two Gaussian functions and the 
Breit-Wigner function convoluted with a Gaussian. The 
Gaussian functions are used to fit the uj and p mesons 
while the Breit-Wigner convoluted with a Gaussian is 
used to fit the p meson. Figure [3] shows the results. 

The parameters of the fit are fixed to the world average 
values of the masses and widths of the three mesons 
and the widths of the Gaussian distributions, which ac¬ 
count for the detector mass resolution, are constrained by 
the values obtained from the simulation. The p meson 
mass resolution is ~85 MeV/c^. Because the invariant 
mass peak of the p meson is partly resolved in the plot, 
while UJ and p meson peaks cannot be resolved, only an 
estimate of the combined yield of these two mesons is al¬ 
lowed. In the fit, the ratio of the w and p mesons, Np/N^^, 
is set to 0.58, derived as the ratio of their corresponding 
production cross sections, <Jpl<Juj = 1.15±0.15, consistent 
with values found in jet fragmentation , multiplied by 
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FIG. 3. (color online) The fitted unlike-sign dimuon spectra 
after background subtraction. 


the ratio of their branching ratios [4l|. 

The yield extraction is performed in bins of px over the 
range 1 < pT < 7 GeV/c for the rapidities 1.2 < \y\ < 
2.2, in bins of y for the px range 1 < pt < 7 GeV/c, and 
in different centrality classes. 

The acceptance and reconstruction efficiency (Ae^ec) 
of the muon spectrometers, including the MuID trigger 
efficiency, is determined by passing pythia 6.421 (de¬ 
fault parameters) [i^ generated (j) mesons through a full 
GEANT [ 4 ^ simulation of the PHENIX detector. The 
PYTHIA simulation output is embedded into real d-|-Au 
data events and then reconstructed in the same manner 
as data. Identical cuts to those used in the data analysis 
are applied to this embedded simulation. 



FIG. 4. (color online) The Aeibc as a function of rapidity 
(x-axis) and px (p-axis) for 0 mesons. 

Figure |4] shows the Asrec as a function of px and ra¬ 
pidity for (j) mesons. The Ae^ec reaches a few percent at 
high-pr and rapidly decreases at low momentum limiting 
this study to pr > 1 GeV/c. 


B. cj} ^ e^e measurement 

The dielectron data set for this analysis was recorded 
in 2008 using a MB trigger that required at least one hit 
in each of the BBC detectors with an ERT (EMCal-RICH 
trigger). The ERT requires a minimum energy deposit 
of 0.6 and 0.8 GeV/c in a tile of 2 x 2 EMCal towers 
matched to a hit in the RICH. The trigger efficiencies 
are described in detail in Ref. [45l| . 

A set of quality assurance cuts is applied to the data to 
select good electron candidates and improve the signal to 
background ratio. Event selection requires the BBC col¬ 
lision z vertex to be reconstructed within ±30 cm of the 
center of the interaction region along the beam direction. 
Charged tracks are reconstructed using the DC and pad 
chambers and requiring px > 0.2 GeV/c. Electrons are 
identified mainly by the RICH detector. Furthermore, 
an electron candidate is required to have a good match 
to a cluster in the EMCAL, and the energy, E, of the 
cluster must satisfy the requirement E/p > 0.5, where p 
is the momentum measured by the DCs. The RICH and 
EMCal combined provide an c/tt rejection factor of order 

lOA 

Electrons and positrons reconstructed in an event are 
combined into pairs. The resulting mass spectrum con¬ 
tains both the signal and an inherent background. The 
uncorrelated part of the background is estimated via an 
event-mixing technique (4^ . which combines electrons 
from different events within the same centrality and ver¬ 
tex class. To extract the (p meson raw yield, the sub¬ 
tracted mass spectrum is fitted by a function consisting 
of several contributions. The (p and uj meson peaks are 
approximated by a Breit-Wigner distribution convolved 
with a Gaussian distribution. Parameters of the Breit- 
Wigner part are set to the global averaged values [4l| 
and the width of the Gaussian part is to account for the 
detector mass resolution. Because the production ratio 
of p meson to w meson is assumed to be 1.15, their ratio 
in the fit is given by the ratio of their branching ratios to 
e“'"e“ in vacuum. The p meson yield is then extracted by 
summing up the bin contents in a ±3cr window, where a 
is extracted from the fit, and subtracting the background 
determined by integrating the polynomial over the same 
window. 

Figure [S] shows an example of the fit to the mixed event 
subtracted mass spectrum. The mixed events were nor¬ 
malized in the mass range between 0.7 to 1.5 GeV/c^. 
The stability of the results was checked by varying the 
normalization region, and the difference was included in 
the systematic errors. The ratio of real to normalized 
mixed events was found to be flat in the region of interest 
for this analysis indicating the validity of the normaliza¬ 
tion. Complete details describing the estimation of the 
systematic uncertainty on the raw yield extraction by 
varying the background normalization, fitting functions, 
range and counting method can be found in Ref. [4^. 

The detector acceptance and reconstruction efficiency, 
as well as trigger efficiency were determined using a 
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FIG. 5. (color online) pr-integrated unlike-sign e"'"e“niass 
spectrum measured by the PHENIX central arms after sub¬ 
tracting the uncorrelated combinatorial background in the 
minimum-bias d-|-Au collisions. The dotted-line shows the 
p meson contribution and the dashed line is the polynomial 
background, which describes the residual correlated back¬ 
ground. The insert shows the raw spectrum before subtrac¬ 
tion, overlaid with the normalized combinatorial background, 
which is estimated using the mixed-event technique. 


GEANT [4^ based simulation of the PHENIX detector 
that tuned the detector response to a set of characteris¬ 
tics (dead and hot channel maps, gains, noise, etc) that 
described the performance of each detector subsystem. 
The (j) —>■ e+e” decays were generated with a realistic 
Pt distribution of (p mesons, within |j/| < 0.5 in rapidity, 
and in full azimuth. The simulated p meson mean and 
width of the invariant mass peak were verified to match 
the measured values in real data at all px- More details 
are given in Ref. [s^ ■ 


C. Calculation of invariant yields and 
nuclear-modification factors 

The dilepton (p meson invariant yield in a given rapid¬ 
ity, Pt, and centrality bin is 

Bii d^N _ 1 1 cN^ 

27rpT dydpT 2 ttpt Ay Apr EtotfVevt ’ 

where Bu is the p —>■ l^l~ branching ratio, iV^ is the 
measured p meson yield, fVevt is the number of sampled 
MB events in the given centrality bin. Ay is the width 
of the rapidity bin, Apx is the width of the pr bin, , 
and c is the BBC bias correction factor [U \^. etot 
includes trigger efficiencies, acceptance and reconstruc¬ 
tion efficiency [i^ . The data points are corrected to 

account for the finite width of the analyzed px bins . 

To gain insight into nuclear medium effects and parti¬ 
cle production mechanisms in d-l-Au collisions, the ratio 
of the p meson yields in d-l-Au collisions to p+p collisions 
scaled by the number of nucleon-nucleon collisions in the 


d-l-Au system, A^coU HH, is calculated as: 

d'^NdAu/dydpx 

where d^NdAn/dydpx is the per-event yield of particle 
production in d-l-Au collisions and d^Npp/dydpx is the 
per-event yield of the same process in p+p collisions. The 
p+p invariant yield used in the RdAu calculation for the 
pA decay channel is the p+p differential cross section 
from Ref. divided by the p+p total cross section, 
42.2 mb. The measured p+p differential cross sections 
for e+e” and K'^K~ decay channels are consistent 
and the p+p reference used in the e“'"e“ RdAu calculations 
is extracted from a Tsallis function fit of the combined 
e“'"e“ and K^K~ spectra. 


D. Systematic uncertainties 

The systematic uncertainties associated with the mea¬ 
sured invariant yields and nuclear modification factors 
can be divided into three categories based upon the ef¬ 
fect each source has on the measured results. All un¬ 
certainties are reported as standard deviations. Type-A: 
point-to-point uncorrelated uncertainties that allow the 
data points to move independently with respect to one 
another and are added in quadrature with statistical un¬ 
certainties. Type-B: point-to-point correlated uncertain¬ 
ties that allow the data points to move coherently within 
the quoted range. Type-C: correlated uncertainties that 
allow the data points to move together by a common 
multiplicative factor, a global uncertainty. 

At forward and backward rapidities, a 5% type-A un¬ 
certainty is assigned to signal extraction, which corre¬ 
sponds to the average variation between the results from 
the different yield extraction fits. Type-B uncertainties, 
point-to-point correlated uncertainties to some degree in 
px, include a 4% uncertainty from MuID tube efficiency 
and 2% from MuTr overall efficiency. An 8% uncertainty 
on the yield is assigned to account for a 2% absolute mo¬ 
mentum scale uncertainty, which was estimated by mea¬ 
suring the J/p mass. 

A 9% (8%) uncertainty is assigned to the —2.2 < y < 
— 1.2 (1.2 < y < 2.2) rapidity bins due to the uncertain¬ 
ties in the Ae-cec determination method itself. The Aerec 
at the lowest px bin is small and sensitive to variations 
in the slope of the input px distribution which affects 
the differential cross section calculations at this px bin. 
To understand this effect, the px dependent cross sec¬ 
tion is fitted by three commonly used fit functions (Hage- 
dorn [^, Kaplan [Elj, and Tsallis [s^) over thepT range, 
2 < px < 7 GeV/c, and the htted functions are extrapo¬ 
lated to the lowest px bin, 1 < pT < 2 GeV jc. The differ¬ 
ences between the values extracted from these fits and the 
measured one at the lowest px bin is within 8%, hence an 
additional 8% systematic uncertainty is assigned to the 
lowest Px bin to account for these differences. For the px 















TABLE I. Systematic uncertainties included in the invariant 
yield and the nuclear modification factor calculations at for¬ 
ward and backward rapidities. As explained in the text, there 


is an 8% type-B systematic uncertainty due to small accep¬ 
tance that impacts the low-pT region only which is not listed 
below. 

Source 

Value (%) 

Type 

Signal extraction 

5 

A 

MuID efficiency 

4 

B 

MuTr efficiency 

2 

B 

Asrec [north/south] 

8/9 

B 

Absolute momentum scale 

8 

B 

Total Type-B 

13 

B 

BBC bias correction for p+p (d+Au) 

10 (0.1-5.8) 

C 

Acoll 

5-7 

C 

TABLE II. Systematic uncertainties included in the 

invari- 

ant yield and the nuclear modification factor calculations at 
midrapidity. 

Source 

Value (%) 

Type 

Signal extraction 

8-15 

A 

ERT trigger efficiency 

1-7 

B 

Acceptance correction 

7 

B 

Electron ID 

9 

B 

Absolute momentum scale 

1-5 

B 

Quadratic sum of (B) 

11-14 

B 

BBC bias correction for p+p (d+Au) 

10 (0.1-5.8) 

C 

Acoll 

5-7 

C 


integrated and rapidity dependent invariant yields the 
8 % uncertainty is assigned to all data bins. Type-B sys¬ 
tematic uncertainties are added in quadrature and are 
shown as boxes. Finally, an overall normalization uncer¬ 
tainty of 10% (0.1%-5.8%) for p+p (d-l-Au) is assigned to 
the BBC bias correction (c) uncertainties and are labeled 
as type-C [2^. The BBC bias correction uncertainties for 
d-l-Au vary with centrality and are considered as type-B 
for centrality dependent spectra. 

For the nuclear modification factor, the type-A sys¬ 
tematic uncertainties arise from the quadratic sum of 
type-A systematic uncertainty in p+p and d-l-Au invari¬ 
ant yields [s^ . Systematic uncertainties including MuID 
and MuTr efficiencies and absolute momentum scale are 
the same between the d-l-Au and p+p invariant yields 
and cancel out. A 9% (7%) systematic uncertainty in 
the Asrec at —2.2 < y < —1.2 (1.2 < y < 2.2), that is 
carried over from p+p, is added in quadrature to type-B 
systematic uncertainties listed in TableHl The acceptance 
limitation at the lowest px bin is the same between d+Au 
and p+p data because it is collected using the same de¬ 
tector and the associated type-B systematic uncertainty 
cancels out. Type-C systematic uncertainties for the nu¬ 
clear modification factor are the quadratic sum of the 


type-C systematic uncertainties for the invariant yields 
of p+p and d-l-Au collisions and the uncertainty associ¬ 
ated with A^coii- The systematic uncertainties are listed 
in TableHl 

At midrapidity, the main contribution to the system¬ 
atic uncertainties of type-A is the uncertainty associated 
with the raw yield extraction, 8%-15%. This uncertainty 
is calculated by varying the background normalization, 
fitting functions, range, and counting methods. Type-B 
systematic uncertainties, point-to-point correlated uncer¬ 
tainties inpr, include uncertainties in the ERT trigger ef¬ 
ficiency (l%-7%) and electron identification (9%), which 
are estimated by varying the analysis cuts. The accep¬ 
tance correction uncertainty is 7% and the momentum 
scale uncertainty is l%-5%. An overall normalization un¬ 
certainty, type-C, was assigned for the p+p (d-|-Au) BBC 
bias correction, which amounts to 10% (0.1%-5.8%). 

For the nuclear-modification factor, the electron iden¬ 
tification and branching ratio uncertainties are the same 
between p+p and d-|-Au and cancel out. As for the rest 
of the systematic uncertainties, the uncertainties for each 
type are added in quadrature between p+p and d-l-Au. 
Additionally, the A"coii uncertainty is added in quadra¬ 
ture to the rest of type-C uncertainties. The systematic 
uncertainties are listed in Table HIl 
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FIG. 6. (color online) 4> meson invariant yields as a function 
of pt in the Au-going direction (solid blue squares) and in the 
d-going direction (solid red points). The cf) meson invariant 
yields at midrapidity are shown as solid green triangles. The 
vertical bars represent the statistical uncertainties and the 
boxes represent type-B systematic uncertainties. The ±6% 
global uncertainty is the associated type-C uncertainty. 


IV. RESULTS 

Figure [6] shows the (j) meson invariant yields in d-|-Au 
and p+p as a function of px in the Au-going direction 
—2.2 < y < —1.2, in the d-going direction 1.2 < y < 2.2, 
and at midrapidity \y\ < 0.35. The invariant yields mea- 
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sured at forward and backward rapidities for all centrali¬ 
ties are significantly different from each other. At midra¬ 
pidity, the magnitude of the invariant yield is close to 
that of the backward rapidity, in contrast to what was 
observed in p+p collisions |38l | where the forward and 
backward rapidities have the same invariant yields and 
both are smaller than the midrapidity invariant yield. 

Figure [7] shows the invariant yields as a function of 
Pj- for the different centrality intervals measured in p+p 
and d-|-Au collisions. The d-l-Au data are presented in 
different centrality classes: 0%-20%, 20%-40%, 40%- 
60%, and 60%-88%, which are determined using the BBC 
over the pseudorapidity range 3.1 < |? 7 | < 3.9, for the 
1.2 < \y\ < 2.2 and \y\ < 0.35 regions. The three panels 
describe the three rapidity regions. The same centrality 
definition is used as in Ref. [^. The wide rapidity cov¬ 
erage of the full PHENIX detector over an extended pr 
range is demonstrated. 

Figure |S] shows the RdAu for </> mesons measured in all 
centralities as a function of px at forward and backward 
rapidities and at midrapidity, also compared with RdAu 
measured in the (j) —>■ K^K~ decay channel [s^. The (j) 
meson RdAu behavior measured in the three different ra¬ 
pidity regions is not the same. In the Au-going direction, 
the RdAu is consistent with unity at the lowest measured 
Pt and increases to RdAu ^ 2 at pr = 2.7 GeV/c. It then 
decreases to unity at the highest measured px ■ In the d- 
going direction strong suppression is observed at the low¬ 
est measured px and then the RdAu increases to unity 
at higher px- At midrapidity, where the measurement 
starts at px =0, the yield of the 4> meson is suppressed 
below 1 GeV/c and RdAu is consistent with unity above 
that. The RdAu measured in the (j) —>■ AT'*'AT “channel 
agrees with the measurement of ^ ^ e“'"e“ within uncer¬ 
tainties in the px region where the measurements over¬ 
lap. The enhanced yield observed in the Au-going di¬ 
rection at intermediate px is characteristic of the Cronin 
effect enhancement of high px particle production 
in d-|-Au collisions relative to scaled p+p collisions with 
a corresponding depletion at low px- The (j) enhance¬ 
ment in the Au-going direction and the suppression in 
the d-going direction are consistent with what is observed 
by ALICE in p-|-Pb collisions at y/s^ =5.02 TeV in 
-4.46 <y< -2.96 and 2.03 <y< 3.53 [H. 

Figure m shows the RdAu of the (j) meson as a function 
of rapidity, summed over the px range 1 < pt < 7 GeV/c 
and integrated over all centralities. The RdAu is en¬ 
hanced in the Au-going direction, shows no modification 
at midrapidity and is suppressed in the d-going direction. 

Figure [TUI shows that the (j) meson nuclear modification 
factor, measured as a function of rapidity, is consistent 
with that of leptons coming from decays of heavy flavor 
particles . This similarity is interesting given that 

heavy flavor quark production is expected to be domi¬ 
nated by hard processes over the accessed px range. In 
contrast, one expects a significant contribution to (j) me¬ 
son production from soft processes, particularly at low 
Px where the yield is dominant. This backward/forward 


enhancement/suppression in d-|-Au collisions is also ob¬ 
served in d-|-Au charged hadron density results measured 
by PHOBOS [l^. The charged hadron result is often 
considered as a rapidity shift in the Au-going direction 
via soft processes. 

Figure [TUI also shows the J/ip RdAu [El]- In the case of 
the J/ip^ the relative RdAu modification as a function of 
rapidity is similar to that of the (j) meson and heavy flavor 
decay leptons. However, the J/ij} suffers from additional 
suppression at backward and midrapidity, which could 
be due to a larger break up cross section or effects in the 
higher-energy-density backward-rapidity region. 

At midrapidity, PHENIX measured the nuclear- 
modification factors of protons and of several mesons 
from the light (tt) meson up to the (j) meson with a 
mass similar to that of protons. The results show no 
significant dependence of RdAu on the mass of the par¬ 
ticle [ 2 ^, whereas the measurements reveal a significant 
dependence of RdAu on the number of valence quarks. 
All mesons show no or very small enhancement in com¬ 
parison to protons in the low to medium px region. 

EigurefTT] shows the RdAu integrated over the px range 
I < Pt < 7 GeV/c as a function of centrality for three 
rapidity intervals: —2.2 < y < —1.2, 1.2 < y < 2.2 
and |y| < 0.35. The RdAu is integrated over the px 
range I < pt < 7 GeV/c and the rapidity ranges 
—2.2 < y < —1.2, 1.2 < y < 2.2 and |p| < 0.35. In the 
Au-going direction, the RdAu shows an indication of en¬ 
hancement with NcoW and in the d-going direction shows 
a strong suppression at high NcoW- At midrapidity, the 
RdAu for 4> mesons reconstructed in e+e" is consistent 
with unity at all centralities. This is also consistent with 
the measurement done in the K~ decay channel. 

Eigure [TUI shows RdAu as a function of px in different 
d-l-Au centrality classes for three rapidity regions cov¬ 
ered by the PHENIX detector: backward, mid, and for¬ 
ward rapidities. In the most peripheral collisions (60%- 
88%), shown in the lowest panel of the figure, the nuclear- 
modification factor measured in all three rapidity regions 
is consistent with unity at all measured px- In the more 
central collisions the RdAu remains at unity only in the 
midrapidity region for px above I GeV/c. This is con¬ 
sistent between the measurements done in the e“'"e“ and 
in the K^K~ decay channels. In the Au-going direction, 
the RdAu reaches a maximum at px close to 2.7 GeV/c 
with an amplitude that grows towards more central col¬ 
lisions. At higher px the RdAu diminishes, approaching 
unity in all measured centrality classes. In the d-going 
direction, all points above 2.7 GeV/c are consistent with 
unity. Below px =2.7 GeV/c Eig. [TUI shows RdAu < I, 
with more suppression in the most central collisions com¬ 
pared to other centralities. 
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FIG. 7. (color online) Invariant pr spectra of the (f> meson 
for different centrality classes in d+Au and p+p collisions 
-^/^NN ~ 200 GeV [ 3 ^ . The vertical bars represent the 
statistical uncertainties and the boxes represent type-B sys¬ 
tematic uncertainties. Type-C systematic uncertainties are 
0.1%-5.8% for d-l-Au invariant yields and 10% for p+p in¬ 
variant yields. The spectra are scaled by arbitrary factors for 
clarity. 



FIG. 8. (color online) The (j) meson nuclear-modification fac¬ 
tor, RdAu, as a function of pr- The solid blue squares indi¬ 
cate the Au-going direction and the solid red circles indicate 
the d-going direction. The upright black triangles are for 
<j) -A- K^K~a.t midrapidity [3^ while the inverted triangles 
are for (j> —>■ e"'"e“. The vertical bars represent the statisti¬ 
cal uncertainties and the boxes represent type-B systematic 
uncertainties. The ±13% global uncertainty is the associated 
type-G uncertainty. 
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FIG. 9. (color online) The (j) meson nuclear modification fac¬ 
tor, RdAu, as a function of rapidity. The solid blue squares 
indicate the Au-going direction while the d-going direction is 
shown in solid red circles. At midrapidity, the upright black 
triangles are for (j) K'^K~ while the inverted trian¬ 

gles are for 0 —>■ e'’'e“. The midrapidity points are slightly 
displaced from zero for clarity. The vertical bars represent 
the statistical uncertainties and the boxes represent type-B 
systematic uncertainties. The ±13% global uncertainty is the 
associated type-C uncertainty. 
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FIG. 10. (color online) JA/r (full green circles), heavy 
flavor decay leptons \2i 1^ (full blue triangles) and (j) me¬ 
son (full red squares) nuclear modification factors, RdAu, as a 
function of rapidity. The vertical bars represent the statistical 
uncertainties and the boxes represent type-B systematic un¬ 
certainties. The type-C systematic uncertainties associated 
with heavy-flavor and J/A meson measurements are 10% and 
8%, respectively. 
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cal uncertainties and the boxes represent type-B systematic 
uncertainties. The ±10% global uncertainty is the associated 
type-C uncertainty. 
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The enhancement at backward rapidity is consistent 
with nuclear px broadening at moderate pr and gluon 
antishadowing, while the suppression at forward rapid¬ 
ity may suggest gluon shadowing and/or partonic energy 
loss. The RdAu enhancement (suppression) at backward 
(forward) rapidity in the most central collisions decreases 
gradually from central to peripheral collisions where for 
the most peripheral collisions RdAu shows no overall mod¬ 
ification. Whether the meson is dominated by soft or 
hard processes remains an open question. 


V. SUMMARY AND CONCLUSIONS 

PHENIX has measured the production of (j) mesons in 
d-l-Au collisions at y/s^ = 200 GeV via their decay to 
dimuons and dielectrons. The (j) meson is measured in 
the forward, d-going direction, 1.2 < y < 2.2 in the pr 
range from 1 to 7 GeV/c, at midrapidity \y\ < 0.35 in 
the pt range below 7 GeV/c, and in the backward, the 
Au-going direction, —2.2 < y < —1.2 in the pt range 
from 1 to 7 GeV/c. 

The measurements reveal that the (j) meson yields in 
d-|-Au compared to binary collision scaled p+p at pT > 
2 GeV/c are suppressed in the d-going direction and en¬ 
hanced in the Au-going direction. The yield measured 
at midrapidity is consistent with binary collision scaled 
p+p. No significant modification of (f> meson production 
is observed in the most peripheral d-|-Au collisions. With 
centrality increasing from semi-peripheral events to cen¬ 
tral the RdAu shows more pronounced suppression in the 
d-going and more increasing enhancement in Au-going 
direction. In that rapidity region, the RdAu has a pr 
dependence, which is characteristic of a Cronin-type ef¬ 
fect. A similar enhancement (suppression) was observed 
by ALICE in p-|-Pb collisions at = 5.02 TeV [^ . 

The rapidity dependence of the RdAu in (/-meson produc¬ 
tion is similar to the open heavy flavor modification 
which may indicate a general rapidity shift compared to 
the symmetric system and/or similar cold-nuclear matter 
effects are present in both the (j) meson and open heavy- 
flavor production. 
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APPENDIX 

Tables IIIIHVll show the numerical values of the mea¬ 
sured invariant yields and RdAu that are plotted in 
Figs. As noted in each caption, the first uncer¬ 

tainty is statistical and the second is type-B systematic. 


ACKNOWLEDGMENTS 

We thank the staff of the Collider-Accelerator and 
Physics Departments at Brookhaven National Labora¬ 
tory and the staff of the other PHENIX participating 
institutions for their vital contributions. We acknowl¬ 
edge support from the Office of Nuclear Physics in the 
Office of Science of the Department of Energy, the Na¬ 
tional Science Foundation, Abilene Christian University 
Research Council, Research Foundation of SUNY, and 
Dean of the College of Arts and Sciences, Vanderbilt Uni¬ 
versity (U.S.A), Ministry of Education, Culture, Sports, 
Science, and Technology and the Japan Society for the 
Promotion of Science (Japan), Conselho Nacional de De- 




13 


TABLE III. Invariant yields of (j> meson production as a function of pr at different d+Au centrality classes. The first uncertainty 
is statistical and the second is type-B systematic. 


min 

max 



Bn dP-N /p YVi-^ 



(GeV/c) 

0%-100% 

0%-20% 

27rp'j-' dydprp ' ’ / ) 

20%-40% 

40%-60% 

60%-88% 

1.0 

2.0 

(4.35±0.41±0.65)xl0“’’ 

(1.01±0.12±0.15)xl0“® 

-2.2 <y < -1.2 
(6.56±1.32±0.98)xl0“’' 

(4.94±0.83±0.74)xl0“’’ 

(1.89±0.32±0.28) x 10“’’ 

2.0 

2.5 

(8.70±0.75±1.13)xl0“® 

(1.95±0.19±0.25)xl0“'^ 

(1.12±0.17±0.15)xl0“'^ 

(5.47±0.87±0.71)xl0“® 

(2.32±0.44±0.30) x 10“® 

2.5 

3.0 

(3.09±0.25±0.40)xl0“® 

(7.35±0.68±0.95) x 10“* 

(3.68±0.62±0.48) x 10“* 

(2.23±0.32±0.29)xl0“® 

(7.05±1.51±0.92)xl0“° 

3.0 

4.5 

(2.70±0.22±0.35)xl0“® 

(5.95±0.82±0.77) x lO”® 

(3.79±0.49±0.49) x lO”® 

(1.87±0.29±0.24)xl0“® 

(6.52±3.22±0.85) x 10“’° 

4.5 

7.0 

(6.35±1.72±0.83) X 10““ 

(1.10±0.39±0.14)xl0“^° 

(9.61±3.51±1.25)xl0““ 

(4.28±2.40±0.56) X 10““ 

(2.59±1.04±0.34) x 10“’’ 

1.0 

2.0 

(2.00±0.25±0.30)xl0“’' 

(3.37±0.65±0.51) x 10“'^ 

1.2 < y < 2.2 
(2.93±0.61±0.44) x 10“’’ 

(3.47±0.76±0.52)xl0“’’ 

(1.82±0.31±0.27)xl0“’’ 

2.0 

2.5 

(3.80±0.54±0.49)xl0“® 

(9.30±1.59±1.21)xl0“* 

(5.47±0.73±0.71)xl0“* 

(4.42±0.71±0.57)xl0“® 

(2.07±0.44±0.27) x 10“° 

2.5 

3.0 

(1.68±0.18±0.22)xl0“® 

(3.47±0.56±0.45) x 10“* 

(2.18±0.37±0.28) x 10“* 

(1.60±0.25±0.21)xl0“® 

(6.52±1.07±0.85)xl0“° 

3.0 

4.5 

(1.48±0.14±0.19)xl0“® 

(2.60±0.63±0.34) x 10“® 

(1.80±0.60±0.23) x 10“® 

(1.44±0.24±0.19)xl0“® 

(6.24±1.33±0.81) x 10“’° 

4.5 

7.0 

(6.23±1.51±0.81)xl0““ 

(1.16±0.27±0.15)xl0“^° 

(7.13±3.73±0.93) x 10““ 

(5.55±1.78±0.72)xl0““ 

(2.39±1.40±0.31) x 10“’’ 

0.00 

0.25 

(9.81±2.88±3.01)xl0“® 

(1.93±0.76±0.55) x 10“® 

\y\ < 0.35 

(1.23±0.53±0.33) x 10“® 

(8.92±4.22±2.42) X 10“® 

(4.87±2.74±1.26) x 10“° 

0.25 

0.50 

(7.18±0.74±1.57)xl0“® 

(1.23±0.22±0.23) x 10“® 

(1.03±0.16±0.22)xl0“® 

(5.46±1.02±1.05)xl0“® 

(4.28±0.81±0.81)xl0“° 

0.50 

0.75 

(3.49±0.35±0.71)xl0“® 

(6.11±1.00±1.30)xl0“® 

(4.52±0.69±0.90) x 10“® 

(3.96±0.54±0.69)xl0“® 

(2.18±0.36±0.39) x 10“° 

0.75 

1.00 

(2.19±0.18±0.42)xl0“® 

(3.54±0.46±0.70) x 10“® 

(3.13±0.36±0.58) x 10“® 

(1.54±0.20±0.25)xl0“® 

(1.28±0.17±0.22) x 10“° 

1.00 

1.25 

(1.18±0.09±0.21)xl0“® 

(2.16±0.24±0.40) x 10“® 

(1.51±0.17±0.27)xl0“® 

(7.44±1.13±1.21)xl0“’’ 

(5.70±0.93±0.93) x 10“’’ 

1.25 

1.50 

(6.20±0.49±1.08)xl0“’' 

(1.26±0.14±0.22) x 10“® 

(8.92±0.99±1.51)xl0“’’ 

(4.87±0.67±0.77)xl0“’’ 

(2.37±0.38±0.39) x 10“’’ 

1.50 

1.75 

(2.75±0.24±0.46)xl0“’' 

(6.32±0.78±1.08) x 10“'^ 

(3.13±0.45±0.51)xl0“’’ 

(1.90±0.30±0.30)xl0“’’ 

(1.65±0.26±0.27) x 10“’’ 

1.75 

2.00 

(1.59±0.15±0.26)xl0“’' 

(3.31±0.43±0.55) x 10“'^ 

(2.34±0.30±0.37) x 10“’’ 

(1.32±0.21±0.21)xl0“’’ 

(4.11±1.03±0.67)xl0“° 

2.00 

2.50 

(9.22±0.62±1.48)xl0“® 

(2.01±0.18±0.32)xl0“'^ 

(1.29±0.12±0.20) x 10“’’ 

(7.15±0.85±1.13)xl0“® 

(3.60±0.56±0.59) x 10“° 

2.50 

3.00 

(3.19±0.24±0.49)xl0“® 

(7.12±0.71±1.10)xl0“® 

(3.54±0.44±0.53) x 10“® 

(2.84±0.39±0.44)xl0“® 

(1.10±0.22±0.18)xl0“° 

3.00 

4.00 

(6.17±0.49±0.94)xl0“® 

(1.15±0.14±0.18)xl0“* 

(8.92±1.10±1.36)xl0“’’ 

(5.52±0.81±0.88)xl0“® 

(1.94±0.44±0.33) x 10“° 

4.00 

5.00 

(9.45±1.25±1.47)xl0“^° 

(2.27±0.37±0.37) x 10“® 

(1.27±0.26±0.20)xl0“’’ 

(8.74±2.06±1.57) X 10“’° 

(1.87±1.10±0.34)xl0“’° 

5.00 

7.00 

(1.02±0.19±0.17)xl0“^° 

(1.46±0.45±0.25) X 10“^° 

(1.98±0.50±0.03) x 10“’” 

(8.27±3.16±1.59)xl0“’’ 

(4.79±2.74±0.90) x 10“’’ 


TABLE IV. RdAu as a function of rapidity of 0 meson summed over the pr range, 1 < pr < 7 GeV/c for 0%-100% centrality. 
The first uncertainty is statistical and the second is type-B systematic. 


1/min 

ymax 

-^dAu 

-2.2 

-2.0 

1.23±0.12±0.20 

-2.0 

-1.7 

1.30±0.13±0.21 

-1.7 

-1.1 

1.34±0.13±0.21 

-0.35 

0.35 

1.02±0.11±0.07 

1.2 

1.7 

0.90±0.11±0.13 

1.7 

2.0 

0.79±0.07±0.11 

2.0 

2.2 

0.44±0.06±0.06 


TABLE V. RdAu vs 
systematic. 

A^coii of (j> meson at 

1 < Pt < 7 GeV/c. 

The first uncertainty is statistical and the second is type-B 

Centrality 

(iVooll) 

-2.2 <y< -1.2 

RdAu 

1.2 <y <2.2 

\y\ < 0.35 

0%-20% 

15.1±1.0 

1.26±0.21±0.21 

0.73±0.12±0.12 

1.03±0.11±0.09 

20%-40% 

10.2±0.7 

1.25±0.20±0.20 

0.65±0.10±0.11 

1.00±0.10±0.09 

40%-60% 

6.6±0.4 

1.18±0.20±0.19 

0.80±0.11±0.13 

0.86±0.09±0.08 

60%-88% 

3.2±0.2 

1.24±0.22±0.20 

0.95±0.17±0.16 

1.04±0.11±0.09 
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TABLE VI. RdAu and pt at different d+An centrality classes. The first uncertainty is statistical and the second is type-B 
systematic. 


(GeV/c) 

0%-100% 

0%-20% 

RdAu 

20%-40% 

40%-60% 

60%-88% 

1.0 

2.0 

0.92±0.10±0.12 

1.07±0.14±0.14 

-2.2 <y< -1.2 

1.03±0.22±0.13 

1.20±0.22±0.16 

0.94±0.17±0.12 

2.0 

2.5 

1.57±0.16±0.19 

1.78±0.20±0.18 

1.51±0.25±0.15 

1.14±0.19±0.11 

1.00±0.20±0.10 

2.5 

3.0 

2.19±0.22±0.26 

2.63±0.29±0.26 

1.95±0.35±0.19 

1.83±0.28±0.18 

1.19±0.26±0.12 

3.0 

4.5 

1.66±0.17±0.20 

1.84±0.28±0.18 

1.74±0.25±0.17 

1.33±0.22±0.13 

0.95±0.47±0.10 

4.5 

7.0 

1.03±0.30±0.12 

0.90±0.33±0.09 

1.17±0.44±0.12 

0.80±0.46±0.08 

1.00±0.41±0.10 

1.0 

2.0 

0.42±0.06±0.05 

0.36±0.07±0.05 

1.2 <y< 2.2 

0.46±0.10±0.06 

0.84±0.19±0.11 

0.91±0.16±0.12 

2.0 

2.5 

0.69±0.10±0.08 

0.85±0.15±0.08 

0.74±0.11±0.07 

0.92±0.16±0.09 

0.89±0.20±0.09 

2.5 

3.0 

1.20±0.14±0.13 

1.24±0.21±0.12 

1.15±0.21±0.12 

1.31±0.22±0.13 

1.10±0.19±0.11 

3.0 

4.5 

0.91±0.10±0.10 

0.80±0.20±0.08 

0.83±0.28±0.08 

1.02±0.18±0.10 

0.91±0.20±0.09 

4.5 

7.0 

1.01±0.26±0.11 

0.95±0.24±0.10 

0.87±0.46±0.09 

1.04±0.35±0.10 

0.92±0.55±0.09 

0.00 

0.25 

0.70±0.29±0.18 

0.69±0.33±0.18 

\y\ < 0.35 
0.65±0.33±0.16 

0.73±0.40±0.19 

0.82±0.51±0.21 

0.25 

0.50 

0.72±0.18±0.08 

0.62±0.17±0.07 

0.76±0.21±0.08 

0.63±0.18±0.07 

1.00±0.28±0.11 

0.50 

0.75 

0.74±0.18±0.07 

0.65±0.19±0.06 

0.71±0.19±0.07 

0.97±0.24±0.10 

1.10±0.29±0.11 

0.75 

1.00 

0.71±0.15±0.07 

0.58±0.14±0.06 

0.75±0.16±0.07 

0.57±0.12±0.06 

0.98±0.21±0.10 

1.00 

1.25 

0.86±0.20±0.09 

0.80±0.20±0.08 

0.82±0.20±0.08 

0.63±0.16±0.06 

0.99±0.26±0.10 

1.25 

1.50 

0.89±0.22±0.09 

0.91±0.23±0.09 

0.95±0.24±0.09 

0.81±0.21±0.08 

0.81±0.22±0.08 

1.50 

1.75 

0.75±0.18±0.07 

0.87±0.22±0.08 

0.63±0.16±0.06 

0.60±0.16±0.06 

1.10±0.28±0.10 

1.75 

2.00 

0.82±0.19±0.08 

0.86±0.22±0.08 

0.89±0.22±0.09 

0.78±0.21±0.08 

0.50±0.17±0.05 

2.00 

2.50 

1.20±0.26±0.12 

1.30±0.30±0.13 

1.20±0.28±0.12 

1.10±0.26±0.10 

1.10±0.29±0.11 

2.50 

3.00 

1.30±0.29±0.13 

1.50±0.34±0.14 

1.10±0.26±0.10 

1.30±0.33±0.13 

1.10±0.31±0.10 

3.00 

4.00 

1.20±0.27±0.13 

1.20±0.28±0.12 

1.30±0.32±0.13 

1.30±0.33±0.13 

0.94±0.30±0.09 

4.00 

5.00 

1.30±0.33±0.14 

1.60±0.43±0.17 

1.30±0.39±0.14 

1.40±0.47±0.15 

0.63±0.40±0.07 

5.00 

7.00 

1.80±0.54±0.21 

1.30±0.52±0.15 

2.70±0.91±0.30 

1.70±0.80±0.19 

2.10±1.30±0.23 
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